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PATENT APPLICATION 



Docket No. : D349 

Inventor (s): Gee L. Lui, Kuang Tsai 

Title: Gaussian Minimum Shift Keying (GMSK) Precoding Communication 
Method 

SPECIFICATION 

Statement\of Government Interest 

The invention was made with^Government support under contract 
No. F04701-93-C-0094 by the Department of the Air Force. The 
Government has certain rights in theVinvention. 

Field of the Invention 

The invention relates to the field of continuous phase 
modulation communications. More particularly, the present 
invention relates to a Gaussian minimum shift keying communication 
method for communicating a preceded digital data stream for 
improved performance. 
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Background of the Invention 



Digital communication systems transmit data by various carrier 
modulation techniques. The spectrum of a digital signal can be 
controlled and made compact by envelope filtering or phase domain 
filtering. Envelope filtering filters the baseband data stream 
prior to upconversion to the carrier frequency and power 
amplification. Envelope filtering for controlling out of band 
power of the transmitted signal, can be found in many operational 
communication systems. However, the power amplifier in the 
transmitter must be linearized or backed off to prevent spectral 
regrowth of the filtered signal. The more efficient phase domain 
filtering approach controls the signal spectrum by frequency 
modulating the filtered signal onto a carrier frequency to form a 
continuous phase modulated (CPM) signal. The CPM signal has a 
constant envelope so that the power amplifier can be operated at 
maximiam output power without affecting the spectrum of the filtered 
signal . 

Gaussian minimum shift keying (GMSK) is a form of continuous 
phase modulation having compact spectral occupancy by choosing a 
suitable bandwidth time product (BT) parameter in a Gaussian 
filter. The constant envelope makes GMSK compatible with nonlinear 
power amplifier operation without the concomitant spectral regrowth 
associated with non-constant envelope signals. These attributes 
render GMSK an attractive modulation scheme in all high throughput 
frequency division multiple access satellite communication systems 
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where only a limited system bandwidth is available with the 
transmitters operating at maximum power output efficiency. 

Data bits are formatted, for example, by non-return to zero 
(NRZ) formatting prior to Gaussian filtering, carrier modulation 
and power amplification. The formatted data is transmitted within 
data symbols of an M-ary alphabet of M possible data symbols. An 
M-ary GMSK signal is defined by the complex envelope described in 

terms of symbol energy E, symbol period T, carrier phase and 
phase pulse g(t) using a modulation index h and the equally 
probable NRZ data symbols belong to an M-ary alphabet. The GMSK 
phase pulse g(t) originates from a frequency response f (t) of the 
Gaussian smoothing filter with a single-sided 3dB bandwidth B, time 
truncated to a time interval of LT, where L is an integer. For 
Gaussian filters with small BT products, the memory length L is 
approximately an integer greater than or equal to 1/BT. The length 
L is the number of elapsed symbol periods for the GMSK signal to 
accrue a full phase change due to a single input symbol and hence 
represents the memory of the GMSK signal. A GMSK signal with 
memory L greater than one is termed a partial response GMSK signal. 
The GMSK signal is communicated to a GMSK receiver subject to 
interference and additive white Gaussian noise (AWGN) . 

An optimum GMSK receiver for an additive white Gaussian noise 
channel demodulates the received signal by coherent demodulation 
into estimated output data stream using a local carrier reference. 
The receiver demodulates by filtering the received signal using a 
bank of Laurent filters that filter the demodulated received signal 
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into a symbol sequence. A Viterbi decoder searches the symbol 
sequence for the most probable transmitted data sequence as an 
estimate of the original NRZ formatted data stream. A typical 
coherent receiver for 2-ary GMSK signal is based on a pulse 
amplitude modulation (PAM) representation of 2-ary CPM signals 
using Laurent matched filters matched to the amplitude modulated 
pulses in the PAM representation, and further employs the Viterbi 
algorithm to optimally demodulate the symbol sequence. For 2-ary 
GMSK transmitters with BT=l/4 and with a channel Bit Error Rate 
(BER) of 0.01, or more, a GMSK receiver consisting of only two 
matched Laurent filters and a 4-state Viterbi algorithm can nearly 
achieve the performance of coherent binary phase shift Keying 
(BPSK) signaling communications. Amplitude modulated pulses have 
also been extended to PAM representation for 4-ary CPM signals. 

In demodulating 2-ary and 4-ary GMSK signals using the Viterbi 
algorithm, a differential decoder is necessary to resolve data bit 
ambiguities while providing a degraded BER with respect to the 
absolute phase demodulation. It is desirable to reduce the BER. 
For noisy channels, the differential decoder yields poor bit error 
rate performance. The Viterbi algorithm typically employs a 
sliding window in the demodulation process where the width of the 
sliding window represents the demodulation memory or delay. The 

surviving state sequence 0^= (Sn^ S^-i/ • • • / S^-w) produced by the 
sliding window Viterbi algorithm at any stage n depends on all the 

demodulated symbols {dj^; 0 < k < n} prior to that stage, where 
Sj^=dQ+di+ . . . +dj^. The term W is the size of the sliding window that 
is dictated by the memory length L of the underlying GMSK signal. 
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This intrinsic data dependency of the survivor state sequences Un 
disadvantageously requires a differential decoder operation in the 
receiver when deciding on the actual demodulated symbol from 
successive survivors of the Viterbi algorithm resulting in a 
differential bit error rate (BER) degradation. These disadvantages 
are solved or reduced by the present invention. 

Summary of the Invention 
An object of the invention is to improve the bit error rate 
(BER) of continuous phase modulation communication systems. 

Another object of the invention is to provide a preceding 
method for continues phase modulation communication systems. 

Yet another object of the invention is to provide a preceding 
method for a continuous phase modulation and coherent demodulation 
communication system to reduce the BER while eliminating the need 
for differential decoders in a receiver. 

Still another object of the invention is to provide a 
preceding method for a Gaussian minimiam shift keying (GMSK) 
continuous phase modulation and coherent demodulation communication 
system to reduce the BER while eliminating the need for 
differential decoders in a receiver. 

The present invention is directed to a data preceding 
algorithm implemented in a modulator of a transmitter to 
substantially improve the resulting BER performance of the 
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continuous phase modulated (CPM) receivers, such as, Gaussian 
minimum shift keying (GMSK) receivers without the use of 
differential decoders while preserving the spectral occupancy of 
the GMSK signals. The preceding algorithm offers performance 
improvement for both 2-ary and 4-ary coherently demodulated GMSK 
signals. The performance improvement afforded by the preceding 
algorithm is up to 2.5 dB for coherent demodulation of 2-ary and 4- 
ary GMSK signals. The preceding algorithm may be implemented using 
a lookup table in the modulator circuitry without altering the 
desired spectral occupancy of the non-precoded GMSK signals. 

Preceding improves the BER performance for coherent 
demodulators of the 2-ary and 4-ary GMSK signals implemented using 
a pulse amplitude modulated signal subject to the Viterbi 
algorithm. The preceding algorithm encodes the source NRZ data 
symbols {dj^} prior to the GMSK modulation so that the cumulative 
phase of the preceded symbols {an) is identical to the absolute 
phase of the source NRZ symbols at every stage of the Viterbi 

algorithm, that is, Tih (ao+ai+ . . . +an) =7chdn, where h denotes 
modulation index. In the Viterbi algorithm, this preceding process 

produces a set of survivor sequences {Un= (S^, S^-i/ • • • / S^-w) } 
satisfying the condition Sk=d]^ at every stage k, thus making the 
differential decoder operation unnecessary. The preceded symbols 
have the same statistics as the source symbols so that the transmit 
spectrum of the GMSK signal is preserved. The BER is reduced using 
the preceding method for GMSK signals with memory of L. The memory 
L of a GMSK signal is related to the bandwidth time product BT of 
the underlying Gaussian filter by L=1/BT, where B is the single- 
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sided 3dB filter bandwidth in hertz. Depending upon the channel 
bit error rate required/ the preceding method will render a signal 
to noise ratio (SNR) improvement of up to 2.5 dB over the same 
modem that demodulates GMSK signals without preceding. The 
proposed preceding method may also be applied to other 2-ary and 4- 
ary CPM signals. The only change required for these cases is to 
modify the pulse amplitude modulation (PAM) filters of the 

receiver. The preceding method encodes the source symbols {dn) 
prior to the GMSK modulation so that the resulting channel symbols 

{ttn} will render an optimal pseudo symbol sequence requiring no 
differential decoding with improved bit error rates. These and 
other advantages will become more apparent from the following 
detailed description of the preferred embodiment. 

Brief Description of the Drawings 

Figure 1 is a block diagram of a preceded coherent Gaussian 
minimum shift keying (GMSK) communication system. 

Figure 2 is bit error rate (BER) performance graph for 
preceded and non-preceded GMSK communication links using a Gaussian 
filter having a bandwidth time product of BT=l/3. 



/// 
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Detailed Description of the Preferred Embodiment 



An embodiment of the invention is described with reference to 
the figures using reference designations as shown in the figures. 
Referring to Figure 1, a baseband representation of an M-ary 
Gaussian minimum shift keying (GMSK) communication system is 
simplistically shown for convenience to have ideal symbol timing 
and carrier phase synchronization. The GMSK transmitter 10 
comprises a data source 11, a data precoder 12, and a GMSK 
modulator 14. The data source 11 continuously generates M-ary NRZ 
data symbols dn chosen from an M-ary alphabet set {±1, ±3, ... ± (M- 
1)}. These source symbols d^ are then preceded by the data precoder 
12 into a preceded symbol sequence that is in turn modulated by 
the GMSK modulator 14. The GMSK modulator 14 includes a Gaussian 
filter 13, a frequency modulator 15, and a frequency converter 16. 
The Gaussian filter 13 is defined by a bandwidth time product (BT) 
that may be, for example, 1/3, where B is the one sided 3dB 
bandwidth in hertz of the Gaussian filter 13 and T is the data 
symbol duration in seconds. For M-ary GMSK signals with h=l/M, 
both main lobe bandwidth and sidelobe amplitude decrease with a 
decreasing BT. The Gaussian filter 13 provides a Gaussian filter 
output G{t) expressed as an accumulative filter sum response of the 
input sequence of preceded symbols a^- 

Git) = Y,a„f{t-nT) 

n 

/// 
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The term f (t) is the well-known GMSK frequency pulse that is a 
function of the BT product and is essentially zero except over a 
time interval of duration LT, where L is an integer representing 
the memory of the Gaussian filter 13. The memory length L is 

greater than or equal to one (L>1) . The frequency modulator 15 
receives and modulates the Gaussian filter output G(t) by a 
predetermined modulation index h that may be, for example, 1/2. In 
general, lowering the modulation index h while keeping the BT 
product constant will reduce spectral occupancy. Preferably, the 
modulation index is set to h=l/M. The frequency modulator 15 
transforms the Gaussian filter output G(t) into a continuous phase 
modulated baseband GMSK signal Z^lt). The signal amplitude of 
V(2E/T) is taken as one. 



The term g(t), which is the integral of the GMSK frequency 
pulse f{t), is the well-known GMSK phase pulse. 

The baseband GMSK signal Zb(t) is then upconverted by the 
converter 16 using a carrier reference 17 and then transmitted over 
a communication channel 18 subject to additive white Gaussian noise 
(AWGN) and potential interference 19. The transmitted GMSK signal, 
along with noise and interference, is received by a corresponding 
GMSK receiver 20 equipped with a frequency converter 21. The 
converter 21 uses a locally generated carrier reference 22 to 
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downconvert the received RF signal into a baseband signal Zj.(t). 
The received baseband signal Zj.(t) is then processed by a trellis 

demodulator 24 to provide data estimates d^ to a data sink 26. The 
trellis demodulator 24 includes a filter bank 28, a sampler 30, and 

a Viterbi decoder 32 implementing a Viterbi algorithm. 

The received baseband signal Zj.(t) is first filtered by the 
filter bank 28 that consists of F matched filters, where F is at 

most Q=2^'^ for 2-ary GMSK signaling, and at most P=3 • [2"'""'^] ^ for 4- 
ary GMSK signaling. The filters in the filter bank 28 are matched 
to the Laurent amplitude-modulated pulses of the transmitted 

baseband GMSK signal Z^lt), and may be implemented as a matched 
filter bank or an integrate-and-dump type filter bank. The filter 

bank 28 provides filtered signals r^(t) for 0 < k < F-1, which are 
sampled by the sampler 30 at every symbol time instants tn=nT to 
produce discrete sample values rj^^^- These sample values are then 

/\ 

processed by the Viterbi decoder 32 to provide the data estimates d^ 

/\ 

to the data sink 26. In order to produce reliable data estimates d 
n/ the processing of the Viterbi decoder 32 must conform to the 
preceding performed by the data precoder 12 on the data symbols d^ 
at the transmitter 10. The number of matched filters used in the 

filter bank 28 also affects the reliability of the data estimates d 
n • 

For an N symbol long 2-ary data sequence {a^; 0 < n < N-1}, the 
baseband GMSK signal Zt,(t) has a Laurent representation. 
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The term Q=2 is the total number of 2-ary amplitude 
modulated pulses {h)^(t)}, and the term {a^^n^ represents the pseudo- 
symbols relating to the 2-ary data sequence {a^} through radix-2 
digits {kj_} for k defined by a summation over the index i with 
equal to zero. 



Each 2-ary amplitude modulated pulse h]^(t) is related to the 
modulation index h and signal memory L through a generalized phase 

pulse c ( t ) . 




Mt) = ]X:lcit + iT-(l-k,)LT) 



sml7ih-nhg(\t\)]/sm(nh). \t\<LT 



c(0 = 



\t\>LT 



0; 



/// 
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The amplitude modulated pulse hy.{t) is time limited to the 
interval [0,DkT] where Dk=mino<i<L{L (2-k^) -i} , for examples, Do=L+l, 
and Di=L-l. The optimal trellis demodulator 24 for an N syit±>ol 
sequence {an; 0 < n < N-1} corrupted by additive white Gaussian 
noise (AWGN) is one that maximizes 2^ correlation metrics for 0 < m 
< 2«-l. 

The term Zj^(t) denotes the baseband signal associated with the 
m-th possible sequence, and the term z^Ct) denotes the AWGN 
corrupted received baseband signal. Expanding Zj^(t) in Laurent 
representation, these correlation metrics can be expressed in terms 

of the pseudo symbols a|^^^. 

The term a|^"^j| denotes the pseudo symbols associated with the m- 
th possible sequence. By supplying the sampled matched filter 

output r^^n to the Viterbi decoder 32, which produces an optimal 
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pseudo symbol sequence {ao,n'' 0 < n < N-1} that maximizes A, the 
best estimate of the transmit sequence [a^} can then be found. The 
trellis demodulator is simplified by retaining the first F matched 

filters {hk(t); 0 < k < F-1} in the A equation, where F is often 
confined to power of two due to the batch nature of the filter 

duration {D^j^ e.g., D2=D3=L-2, D4=D5=D6=D7=L-3 . For the case of 
F=2, the first two 2-ary matched- filters and the corresponding A 
equation are explicitly given by matched filter equations hg and h^. 

1=1 



c(/-(z-i)r) 



The optimal pseudo symbol sequence produced by the Viterbi 
decoder 32 at any stage n inevitably involves all the demodulated 
symbols prior to that stage. 



a. =J"-a, 



O.n-2 
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The term J=exp{j7Ch) depends on the modulation index h, and is 
J=j for the case of h=l/2. This intrinsic data dependency requires 
a differential decode operation when demodulating the actual data 

symbol a^. 

a„ =-7-«o.«-«o.„-i 

This results in a differential BER degradation comparable to that 
of DPSK as compared to BPSK. 

The purpose of the data precoder 12 is to encode the source 
symbols {dn) prior to the GMSK modulator 13 at the transmitter so 
that the resulting preceded channel symbols {a^^} will render an 
optimal pseudo symbol sequence {ao,n} requiring no differential 
decoding in the receiver, thereby improving the data demodulation 
performance. In mathematical terms, the data precoder 12 carries 
out a data mapping. 

The data mapping must provide a resulting expression for the pseudo 

symbol ao,n involving only d^ and not {di; 0 < i < n} . There is no 
known systematic routine that can be used to determine such a 
mapping. Through repeated trial and error, and essentially by 
chance, one such mapping for 2-ary GMSK signal with modulation 
index h=l/2 has been found. 
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a« = ^„ • d„., = [d„ - d^, + iL^, {a, = d, ) 



This data mapping preserves the transmitted spectrum of the 
GMSK signal because the preceded symbols {a^^} are still equally 
probable as the source symbols {dn}. The data mapping can be 
implemented in the data precoder 12 through the first 2-ary 
precoder lookup table. 



First 2-ary Precoder Lookup Table 


dn 


dn-l 


an=fdn-dn-i+l]niod4 


-1 


-1 


+1 


-1 


+1 


-1 


+1 


-1 


-1 


+1 


+1 


+1 



The data mapping defined by the first 2-ary precoder lookup 
table results in an optimal pseudo symbol sequence produced by the 
Viterbi decoder 32. 

a,,„=r-r'=r'-d„ 

a,„=r-'-J'''-'''-'''-'=r-d„d„.,d„., 
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With the decoding states defined as Sn= (dn, d^-i) / a 2 -state 
2"^-branch Viterbi algorithm is sufficient for demodulating source 
symbols {dn} when two 2-ary matched-filters are used in the filter 
bank 28. 

A second 2-ary data precoding mapping for 2-ary GMSK signal 
with a modulation index h=l/2 has also been found. 

«„=(-ir-^„-<-,=(-i)" k-<-.+iL4 («o=^o) 

The second data precoding mapping also preserves the transmit 
spectrum of the GMSK signal, and can be implemented in the data 
precoder 12 through the second 2-ary precoder lookup table. 



^ Second 2-ary Precoder Lookup Table 




dn-1 


an=[dn-dn-i+l]mod4 
n: even 


an=[dn-dn-i+l]mocl4 
n: odd 


-1 


-1 


+1 


-1 


-1 


+1 


-1 


+1 


+1 


-1 


-1 


+1 


+1 


+1 


+1 


-1 



The data mapping defined by the second 2-ary precoder lookup 
table results in an optimal pseudo symbol sequence produced by the 
Viterbi decoder 32. 
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j d„; n . even 

l.n "1 

d„ ; n : odd 



a, = < 



-j d„ d„_,-d„_^; n. odd 



2 

With the decoding states still defined as 3^= (d^/ d^-i) / a 2 - 
state 2'^-branch Viterbi algorithm is again sufficient for 
demodulating source symbols {dn} when two 2-ary matched filters are 
used in the filter bank 28. The choice between the first and the 
second 2-ary data precoding mapping is arbitrary. Extensive 
simulations of the first and second mappings consistently yield 
identical demodulation performance. 



Extending the 2-ary trellis demodulator to a 4-ary trellis 
demodulator is based on expressing every 4-ary symbol a^eill, ±3} in 
terms of two 2-ary symbols a^^^ and a^^' . 
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This 4-ary an expression enables any 4-ary GMSK signal to be 
expressed as a product of two 2-ary GMSK signals with h^^^=h and 
h^"^^==2h as the modulation indices, respectively. 



Expressing each 2-ary signal constituent into Laurent 
representation and combining the product, the baseband 4-ary GMSK 

signal, for an N symbol long 4-ary data sequence {a^; 0 < n < N-1}, 
takes an amplitude modulation pulse form. The signal amplitude of 

V(2E/T) is taken as one. 



The term P= 3Q is the total number of 4-ary amplitude 
modulated pulses {fk(t)}, and the terms {b^^n^ the 4-ary pseudo 

symbols associated with the 4-ary data sequence {an}. All the 4-ary 
entities can be obtained from respective 2-ary counterparts by 
following a systematic enumeration approach. The optimal 4-ary 
trellis demodulator 24 is identical to that for 2-ary GMSK with the 
following replacements for a A equation and rj^^n equation. 




^ p-\ N-\ 
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The duration of the 4-ary amplitude modulated pulses {f^lt)} is 
also presented in batches, that is, Do=L+l, Di=D2=L, 03=. . .=Di]^=L-l . 
The 4-ary trellis demodulator is simplified by retaining the first 
F=l, F=3 or F=12 matched filters in the filter bank 28. For the 

case of F=3, the 4-ary matched filters and the corresponding A 
equation and pseudo symbols can be explicitly expressed. 

/2 (0 = h^\t) ■ h^\t + 1) = Kit, h = \{t^\,h = /i^')) 
A(w) = Re[ro „ • bl„ + „ • 6*„ + „ • bl„\ 

Kn 
K 



The term J= ( jTih) = (1+j ) /V2 for the case of h=l/4. The (0) and (1) 
terms indicate the modulation index being used for the 2-ary 

constituents, that is, h^°^=h or h^-'-^=2h. As in the 2-ary GMSK 
case, the 4-ary pseudo-symbol bo,n any stage n also involves 



~ "0,n "0,n ~ 

-/*(°) .^0) _ /2oi"+Oo+Ol+-+<»»-i 
~ "0,n-l "0,/i ~ 

,(0) ^(1) _ 7«i'"+Oo+«i+-+a«-i 



"0,n "O.n-1 
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prior demodulated symbols, and must resort to the same differential 
decode operation when demodulating the actual data symbol. Again, 
through trial and error, two spectrum preserving data preceding 
mappings have been found for the 4-ary GMSK signal with modulation 
index h=l/4. 



These preceding mappings can be implemented in the data 
precoder 12 through the 4-ary precoder lookup table. 



/// 
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4-ary Precoder Lookup Table 






' A 


'^n '•^n ^n-l ■'•-'itioao 


-3 


-3 


+ 1 


+3 


-3 


-1 


-1 


+1 


-3 


4-1 


-3 


-1 


-3 


+3 


4-3 


-3 


-1 


-3 


+3 


-3 


-1 


-1 


-fl 


+3 


-1 


+1 


-1 


4-1 


_1 


+3 


-3 


-1 


+1 


_3 


-3 


-1 


+1 


-1 


+3 


-3 


+1 


H-l 


+1 


+3 


+1 


+3 


-1 


+1 


+3 


-3 


-1 


+1 


+3 


-1 


-3 


-1 


+3 


+1 


+3 


-3 


+3 


+3 


+1 


+3 



The first preceding mapping otn= [dn-dn-i+llmodS results in an 
optimal pseudo symbol sequences produced by the Viterbi decoder 32. 



Kn 


=r 




Kn 






Kn 


= j" 
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The term j" belongs to the set {±1, ±j, (±l±j ) /V2 } , and both 
a^^' and a^^^ are deterministic functions of dn and d^-i- Similarly, 
the second precoding mapping C<^n~'^n-1'''3] mode results in an optimal 
pseudo-symbol sequence produced by the Viterbi decoder 32. 



Kn=J 



0) 



(0) 



In both cases, with the decoding state defined as Sn=(dn)/ a 
4''"-state 4^-branch Viterbi algorithm is sufficient for demodulating 
the source symbols {dn} when three 4-ary matched-filters are used in 
the filter bank 28. 

Figure 2 quantified the performance improvement achieved 
through the data precoding for both the 2-ary and 4-ary GMSK 
signals with BT=l/3. Simulation data show that, depending on the 
channel bit-error-rate of interest, a GMSK modem employing data 
precoding will render a 0.5 dB to 2*5 dB signal-to-noise ratio 
(SNR) enhancement over the same modem that employs no data 
precoding. 

The function of the precoder 12 in the GMSK transmitter 10 is 

to precondition the symbol sequence aj^ as an effective reverse 
function of differential encoding that intrinsically results from 
the GMSK modulation process. The precoding produces absolute phase 
demodulation achieved within the GMSK receiver 20. This absolute 
phase demodulation eliminates the need for differential decoding of 

-22- 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 




matched filters 28 while providing an improvement in signal 
detection performance. The preferred precoding algorithms are 
specific to M-ary CPM signals. 

The present invention is directed to the precoding of a data 
sequence into an encoded sequence of transmitted symbols, to avoid 
differential decoding and for improving the BER using Laurent 
filtering. In the preferred form, a precoder is applied to 2-ary 
and 4-ary symbol sets used in a GMSK transmitter having a Gaussian 
filter defines by respective BT products and frequency modulator 
modulation indices. The preferred GMSK receivers included matched 
filters 28, sampling 30, and Viterbi decoding 32. The general form 
of the invention is a precoding method applicable to any M-ary 
symbol set, BT product, modulation index, bank of match filters, 
and Viterbi decoding algorithms. Those skilled in the art can make 
enhancements, improvements and modifications to the invention, and 
these enhancements, improvements and modifications may nonetheless 
fall within the spirit and scope of the following claims. 
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